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ABSTRACT 

A successful attempt was made to analyse about 6000 single pulses of PSR B1133+16 ob- 
tained with the 100-meter Effelsberg radio-telescope. The high resolution (60 micro-seconds) 
data were taken at a frequency of 8.35 GHz with a bandwidth of 1 . 1 GHz. In order to examine 
the pulse-to-pulse intensity modulations, we performed both the longitude- and the harmonic- 
resolved fluctuation spectral analysis. We identified the low frequency feature associated with 
an amplitude modulation at /4 ^ 0.033Pj^', which can be interpreted as the circulation time 
P4 ^ 30f 1 of the underlying subbeam carousel model. Despite an erratic nature of this pulsar, 
we also found an evidence of periodic pseudo-nulls with P4 = 28.44Pi. This is exactly the 
value at which Herfindal & Rankin found periodic pseudo-nulls in their 327 MHz data. We 
thus believe that this is the actual carousel circulation time in PSR Bl 133+16, particularly 
during orderly circulation. 
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1 INTRODUCTION 

Generally pulsars are weak radio sources, and even with the sensi- 
tivity of current radio telescopes single pulses are observable only 
from the strongest objects. Due to the steep spectrum of pulsar 
radio-emission, single pulse observations were readily made at low 
frequencies, where the sensitivity requirements for antenna and re- 
ceiver systems were not an issue. PSR B1133+16 is one of the 
brightest pulsars in the norther n hemisphere (32 mjy at 1.4 GHz, 
spectral index of -1.5 iMaron et al. 2000). Individual pulses of 
PSR B 1133+ 16 were recorded during numerous observing cam- 
paigns, and analyses of single pulse behaviour were performed 
at the observing frequencies from 110 MHz to 1.4 GHz (Backer 
I973I: IXavl or. Manchester & Huguenin 1975; Now akowskili 19961 : 



Weltevrede et al..,2006., 2007. : Herfindal & Rankin.2007h . With the 



recent development of high sensitivity receivers, and given the fact 
that PSR Bl 133+16 is a relatively strong source, measurements of 
its single pulses with high time resolution were possible at frequen- 
cies as high as 8.35 GHz. 

The most common feature revealed by single-p ulse sequences 
is rel a ted to the phenomenon of drifting subpulse s jPrake & CraftI 
1 19681 : ISuttonetalJll970l : 'Weltevrede etalioof). The parameters 
of the drift we r e first desc ribed by Bac ker ( 1973), later modified by 
iGil & Sendv^ i2003h and lGupta et al] l2004h . The first of the four 
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main periodicities involved is the basic pulsar (spin) period P[, i.e. 
the interval between our successive sampling of the pulsar beam. 
The overall pulsar beam is organized into a number of subbeams, 
only some of which cross our sightline during any given pulse, and 
thus directly correspond to the number and the locations of sub- 
pulses we observe. The second periodicity, P2, is the angular dis- 
tance between adjacent subpulses, or the corresponding subbeams, 
as viewed along rotational longitude. The magnetic azimuthal sep- 
aration, say T], between adjac ent subbeams implied by P 2 can be 
estimated (e.g., see Eq. 3 of 'Peshpande & Rankinll200lh . For N 
subbeams, when evenly placed within the emission cone, t] would 
be = 360° /N. The next period P3 (usually measured in units of P[, 
and as per its initial definition) is the distance between adjacent drift 
bands, or in other words, the interval after which the subsequent 
subbeam reaches the same phase as the preceding one. In case of 
fast drifts, however, one has to take into account possible aliasing 
effect, which can cause a large degree of ambiguity. Fluctuation 
spectral analysis can help to resolve the aliasing issue, and is there- 
fore commonly used. In such Fourier analysis, the periodic phase 
modulation associated with the subpulse drift manifests as a spec- 
tral feature at the corresponding fluctuation frequency, /j = I/P3. 
This feature is usually seen at relatively high fluctuation frequen- 
cies, and hence, often called the high frequency feature (HFF). Fi- 
nally, we can introduce the total circulation time of the pattern, or 
the inte rval after which a given subbeam reaches the same phase, P4 
(P3 in iRuderman & Sutherlan j|T975L RS75 henceforth, who in- 
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traduced this feature theoretically). The circulation time becomes 
measurable if the subbeams vary in intensity, as they most often do, 
and such periodic amplitude modulation would be apparent as an 
accompanying spectral feature at frequency = I/P4. This spec- 
tral feature, when present, is invariably located in the lower band 
of the fluctuation spectrum, and is referred to as the low frequency 
feature (LFF). According to the carousel model , first introd uced 
by RS75 and developed bv|Deshpande & RankinI ( Il999ll200ll) and 




iGil & Sendv3 ( I2OO0I l2003h . the two periodicities P3 and P4 are 
connected to each other via a direct relationship P4 = NP3 . 

PSR Bl 133+16 apparently shows the phenomenon of drifting 
subpulses, although it was shown that the organized character of 
the subplulses appears only for some finite timespans, outside of 
which the behaviour of the individual pulses is highly chaotic. Re- 
centlylWeltevrede et al. ( 2007) detected the high-frequency (phase- 
modulation) spectral feature corresponding to period of ^3 = (3 ± 
2)Pi. Moreover, they also detected a low-frequency (amplitude- 
modulation) feature corresponding to period P4 = (30 ± ^)Pi. The 
very large uncertainty in the former estimate follows from a fact 
that the a uthors were not able t o reso lve the issue of aliasing in their 
analysis. iHerfindal & RankirJ ilOOt HR07 henceforth) attempted 
to resolve the aliasing problem and demonstrated that the actual 
Pi = (1.237 ± 0.01 l)Fi. Their spectral and profile-folding anal- 
yses confirmed the existence of the low-frequency feature (LFF) 
in the spectrum, corresponding to a period of P4 = 28.44f 1. The 
unprecedented accuracy of the latter estimate is a result of sophis- 
ticated analysis, including the cart ography technique developed by 
iDeshpande & RankirJ fl999l . l200lh . as well as detection of periodic 
pseudo-nulls of very short duration. It is worth noting that the value 
of N=28.44/1.237 is exactly 23 in their analysis, as it supposed to 
be for a parameter describing the number of subbeams contribut- 
ing to the observed drifting subpulse phenomenon. Thus, the value 
of P4 = 28.44/'i represents the carousel circulation time in PSR 
Bl 133+16 (HR07). 

In this paper we attempt to find the pulse modulation features 
using spectral analysis of our high frequency (8.35 GHz), high time 
resolution (60 /is) single pulse data of PSR B1133+16. 



2 OBSERVATIONS AND DATA ANALYSIS 

Our observations of single pulses of PSR B 1133+ 16 were made 
with the 100-meter radio telescope of Max-Planck Institute for Ra- 
dioastronomy in 2004. We used a cooled secondary-focus receiver 
equipped with HEMT amplifiers, with an observing frequency of 
8.35 GHz and 1.1 GHz bandwidth. The pulsar sign al was recorde d 
by EfFelsberg Pulsar Observations System (EPOS. JessnerllT996h . 
with a time resolution of 60 micro-seconds. To achieve the best 
pulse coverage, we used a gating system, in which 1024 bins were 
used to cover only about 5% of the pulse period around the main 
pulse. Data were gathered during two separate 60-minute sessions 
(c.a. 3000 pulses each). 

Figure 1 shows a sample of the data we gathered; a set of 
256 consecutive individual pulses, with its average pulse profile at 
the top. One can see that only some of the individual pulses were 
clearly detected. For a majority of the pulses, the radio emission 
was very weak or even missing. Naturally, the typical phenomena 
related to subpulse drifting, like the appearance of drift bands or 
high-frequency features in the modulation spectra, may not be read- 
ily apparent. However, contrary to our low expectations, we man- 
aged to detect all periodicities appearing in the pulse intensity mod- 
ulations visible in better quality low-frequency data. The level of 



i 752 



712 



652 



552 



512 



-5 



Longitude (deg) 



Figure 1. A sequence of pulses from 512th to 768th pulse from observation 
B, showing clearly a periodic intensity modulation (conf. Tables 1 and 2). 
The thick line at the top represents the integrated profile for this sequence. 



confidence of these detections varied between different segments of 
the data set. For example, a low-frequency amplitude-modulation 
feature was very prominent in the analysis of one set of pulses, 
while in some sets the feature can be recognized in the spectrum, 
but not at a level sufficient to claim a clear detection. However, the 
particular set (segment B; pulses 5 12-768) shown in Figure 1 is one 
of our best-quality intervals, where the periodic intensity modula- 
tion, with P4 close to 28 Pi, is clearly visible. 



2.1 Longitude-Resolved Fluctuation Spectra 

In order to analyse the gathered data, we used the me thod of 
Longitude-Resolved Fluctuation Spectra (LRFS hereafter, [Backed 
which detects the regularities in the intensity variations 
of individual pulses. To implement the method, we divided the 
pulse series into blocks of 256 pulses, and performed the one- 
dimensional fast Fourier transform (EFT) for each longitude bin, 
and computed the power spectrum. The spectra were then normal- 
ized for each individual longitude bin, resulting in 1024 indepen- 
dent power spectra. These were then summed to form an integrated 
power spectrum separately for each of the 256-pulse stacks. 

A feature with periodicity around 2.85 Pi was identified in 
the fluctuation spectra at all longitudes, as well as in the integrated 
spectrum. This feature was present in all our pulse sequences, and 
it's origin was traced down to the AC power line. We confirmed that 
this feature was not associated with the pulsar signal by obtaining 
LRFS for the off-pulse region. The feature was extremely strong, 
potentially masking any real features present in the contaminated 
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part of the spectram. Therefore, we decided to remove it from the 
fluctuation spectra. To achieve this, we reconstructed a continuous 
time series from the available gated data and calculated its FFT. In 
this fluctuation spectrum, we replaced the values at the locations 
corresponding to the fundamental frequency of the contaminant, 
and its harmonics (separated by the pulsar frequency as a conse- 
quence of the gating), with zeros. An inverse Fourier transform of 
thus edited spectrum yielded a modified time series, free of this in- 
terference, which after the re-application of gating allowed us to 
repeat the LRFS analysis. 



2.2 Harmonic resolved fluctuation spectra 

The LRFS gives us a wealth of information about the fluctuations 
present in the sequence, but also has a serious drawback of aliasing 
due to poor sampling of fluctuations at a given longitude only once 
every pulse period. In order to deal with the issue of aliasing, we 
need to overcome the basic limitation of this sampling. Noting that 
the fluctuations are sampled more than once in the finite duration 
of the pulse, one can combine the fluctuation spectra of the dif- 
ferent longitude bins with appropriate phases corresponding to the 
longitude separation to obtain the Harmonic-Resolved Fluctuation 
Sp ectra (HRFS hereafter ) . This appro ach was used and developed 
bv lDeshpande & RankiiJ ( Il999ll200ll) . In practice, this is achieved 
by Fourier transforming the entire continuous sequence, which can 
be reconstructed (from the gated version) by using the available on- 
pulse samples, and filling the rest of the longitude regions with ze- 
ros. In this process, the pulse sequence is efl'ectively multiplied by 
a periodic window function which leads to convolution of the spec- 
trum associated with the un-gated pulse sequence with the Fourier 
transform of the periodic window (gating) function used. 



3 RESULTS 

We analysed several 256-pulse (c.a. 300 second) data blocks from 
both of our 60-minute observations (we call them A and B for the 
purposes of referencing here). For some of these blocks, we have 
identified a prominent low frequency feature (LFF) in the integrated 
spectra obtained with LFRS, as well as HRFS, analysis. 

Figure 2 shows the result of our LRFS analysis for an 
interference-cleaned 256-pulse stack, and Table 1 lists the values of 
the measured frequency of the LFF and the corresponding peri- 
odicity P4 = I//4 for those data blocks in which these features were 
detected using the LRFS analysis. For the remaining data blocks the 
LFF was either not present at all (no identifiable peaks close to this 
frequency in the power spectrum), or its signal-to-noise ratio was 
too small to claim its presence beyond doubt (usually due to other 
features that appeared in the power spectrum). The values cor- 
responding to the LFF vary between the blocks, but stay consistent 
within the estimated uncertainties. 

For the observation A, it was possible to identify the LFF in the 
average spectrum, obtained as a sum of the spectra from all the data 
blocks. The value of the LFF frequency from the average spectrum 
(see Table 1) was close to the 28.44Pi obtained by HR07 as the 
putative circulation time of their carousel solution. The same was 
not possible for observation S, as the LFF frequency varies signif- 
icantly from block to block, resulting in a smearing of the feature 
in the averaged spectrum. However, it is worth noting that in the 
best block B (512-768, see Fig. 1) with the lowest estimated errors, 
the P4 = (28.4 ± \.2)Pi was also very close to the HR07 value. 



Table 1. Parameters of the observed low-frequency feature (LFF) obtained 
from the LRFS analysis. Table includes also the value obtained from the av- 
erage spectrum of the whole observation A, as well as the weighted average 
of the LFF frequency, and the value of based on it (in parentheses). 



Observation 


Pulse 




LFF 




sequence 


Period 


Freq. [Pj'] 


A 


1-256 


28.4±3.3 


0.035 ± 0.004 


A 


512-768 


28.4±3.5 


0.035 ± 0.004 


A 


2109-2365 


32.0±2.2 


0.031 ±0.002 


B 


1-256 


36.5±4.5 


0.027 ± 0.003 


B 


512-768 


28.4±1.2 


0.035 ± 0.002 


B 


2304-2560 


34.6±9.5 


0.029 ± 0.007 


Freq. weighted average 


(30.9±1.1) 


0.0323 ± 0.0012 


Average spectrum of A* 


28.3±7.2 


0.035 ± 0.009 


* - see text in 


Section|3]for 


details. 





Table 2. Parameters of the observed low-frequency feature (LFF) obtained 
from the HRFS analysis. Table includes also the value obtained from the av- 
erage spectnim of the whole observation A, as well as the weighted average 
of the LFF frequency, and the value of P4 based on it (in paientheses). 



Observation Pulse 




LFF 


sequence 


Period [Pi] 


Freq. [P^'] 


A 1-256 


28.4±3.1 


0.035±0.004 


A 512-768 


32.0±4.6 


0.031±0.004 


A 2109-2365 


32.0±2.3 


0.031±0.002 


B 1-256 


36.5±6.3 


0.027±0.005 


B 512-768 


28.4±3.0 


0.035±0.004 


B 2304-2560 


32.0±8.8 


0.031±0.010 


Freq. weighted average 


(31.5±1.4) 


0.0317±0.0014 


Average spectrum of A* 


25.6±3.3 


0.039±0.005 



* - see text in Section|3]for details. 



Table 1 also contains the weighted average of the LFF period, ob- 
tained from all individual block measurements. 

The result of our HRFS analysis of the same 256-pulse se- 
quence (as used for the LFRS in Fig. 2) is shown in Fig. 3. The 
HRFS for these pulse sequences also show the LFF, usually with 
another feature (at the frequency equal to 1.0 minus the LFF fre- 
quency; conf. Figure 3), which would be expected to be present 
as the LFF's symmetric counterpart (as in a pair of sidebands about 
each of the pulsar harmonics) in case of amplitude modulation. The 
results of the measurement of the LFF parameters (see Table 2) 
agree with those from the LRFS analysis. Similar to the case with 
LRFS, we were able to find the LFF in the average HRFS only for 
the observation A, and the estimated weighted-average frequency 
(and the corresponding period) associated with the LFF are also 
presented in Table 2. 



4 DISCUSSION AND CONCLUSION 

The average fluctuation frequency of the LFF , that we identi- 
fied (in selected segments of our data), corresponds to a period of 
(31.5± 1.4)/'i (see Table 2). We have chosen the value calculated by 
the means of HRFS analysis as the best representative for a num- 
ber of reasons (although one has to note that values from LRFS 
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Figure 2. Longitude-resolved fluctuation spectrum for PSR B1133+16 
computed for 8.35 GHz observations (using the sequence from 1st to 256th 
pulse from observation A). The main panel presents, in grey-scale, the full 
fluctuation spectrum for all longitude bins corresponding to the integrated 
pulse profile, which is shown in the left hand panel. The bottom panel 
presents the integrated fluctuation spectrum; note the prominent feature 
close to 0.03Pj'. This feature is clearly visible also in the LRFS corre- 
sponding to the leading bright pulse component. The white band visible in 
the LRFS greyplot is a result of removal of a whole range of frequencies 
that were distorted by AC power line interference. 



and HRFS agree well within the error estimates). First of them is 
that HRFS adds the contribution associated with the modulations 
potentially coherently, using the implicit connection between mod- 
ulation seen at different longitudes, as against the incoherent addi- 
tion of modulations contribution across longitudes in LRFS. The 
former therefore is expected to provide a higher signal-to-noise ra- 
tio for coherent modulations, than the latter case. Also, while in 
case of HRFS the spectral power is divided between the feature 
and its alias, for amplitude modulated features the same goes for 
LRFS analysis, we just do not look at the negative part of the spec- 
trum. And finally, having two separated, but related features (alias 
frequency is 1 minus LFF frequency) allows for more precise mea- 
surement of the LFF position in the spectrum. 

The value we obtained is consistent with other single-pulse 
studies of PSR Bl 133+16 (see Table 3). More importantly, a few 
segments of the best quality data show P4 close to 28.44Pi, the 
value given by HR07. They estimated this refined value of us- 
ing profile folding over the circulation period, which is desirably 
sensitive to periodic pseudo-nulls. 

Our data do not show a direct evidence of the HFF correspond- 
ing to subpulse drift because of the erratic nature of the individual 
pulses (see Fig. 1); this is quite normal behaviour for this pulsar, 
which is - at our observing frequency of 8.35 GHz - amplified by 
the fact, that most of the incoming pulses remain below our detec- 
tion threshold. Since P3 (aliased or not) measures the distance be- 
tween subpulse drift bands, it is clear that in our data, where most 
of the pulses are missing, this periodicity cannot be easily detected. 
In addition one has to remember, that even at lower observing fre- 
quencies, the drifting subpulse phenomenon is not always present. 
However, one should not ignore a putative drift feature at frequency 
about 0.81/'7',that we found in selected segments of our data (hid- 
den in a forest of features in the average HRFS spectrum), which 
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Figure 3. Harmonic-resolved fluctuation spectrum for PSR B1133+16 
computed for the 8.35-GHz observations. The grey-scale plot in the central 
panel shows the full spectrum, and the bottom panel - the integrated spec- 
trum. In both of them, two features clearly stand out - one close to 0.03P|"' , 
which is the low frequency feature (LFF), and its symmetric counterpart at 
~ 0.97/"!^'. The contributions, shown in the main panel, corresponding to 
any given frequency appear stretched along the harmonic-number axis. This 
is an anticipated consequence of periodic windowing, including the natural 
pulse-shape mask and any deliberate windowing, such as gating. Left panel 
shows the relative strengths of pulsar frequency harmonics (equivalent to 
the Fourier transform of the average pulse profile). The presented spectrum 
was obtained for the same sequence of 256 pulses, as used for the LFRS in 
Fig. 2 (i.e. the first 256 pulses from observation A). One has to remember 
that, for the purpose of HRFS analysis, we had to reconstruct a continuous 
time series from the available gated observations. 



corresponds to the periodicity of about = 1.237/'i. We note 
its similarity to Fig. Al in HR07, where this drift feature at fre- 
quency about 0.81 P^' is clearly present. It is also worth noting 
that /'4/P3 = 28.44/(1.23 ± 0.01) is exactly 23, which represents 
the number of circulating beamlets in the carousel model of HR07. 

There is no fundamental reason why the LFF corresponding 
to amplitude modulation also should be weak or even missing in 
our data. Even when the sub-pulse phase behavior is less ordered, 
and the primary phase fluctuation feature (HFF) may not be ap- 
parent in the spectrum, the low-frequency feature indicative of cir- 
culation period may stand out clearly representing the overall in- 
tensity variation across magn etic azimuth. The first of such exam- 
ples of LFF was discussed bv lAsgekar & Deshpandd ( l200lLl2005l) 
based on their low-frequency (35 MHz) data of B0943+10 and 
B0834+06. Also, a very good example is the so-called Q-mo de 
in PSR B0943+10 analysed bv lRankin & Sulevmanov3 ( I2OO6I) in 
their 327-MHz Arecibo data. Unlike in the regular B-mode of this 
pulsar with driftbands clearly visible and the HFF appearing in the 
fluctuation spectrum, single pulses in the Q-mode are completely 
chaotic, and no driftbands app ear in pulse sequences (see Fig. 4 in 
iRankin & S ulevmanovall200 (;). Still, in this mode, the LFF corre- 
sponding to carousel circulatio n time is clearly visible in th e fluc- 
tuation spectrum (see Fig. 6 in lRankin & Sulevmanovall2006 ). Yet 
another example of pulsars with erratic subpulses, but nevertheless 
showing the LFF that can be interpreted as the carousel circulation 
time, is PSR B0656+14 (see footnote^. 

It is interesting and important that we are able to identify the 
LFF in our erratic data using both, LRFS as well as HRFS analysis. 
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Table 3. Results of observed Low Frequency Feature in PSR Bl 133+16 
found in the literature and from our analysis 



References 


Observation 
frequency 


Period from LFF 
(in units of Pi) 


Clan jp,- AiTiiiiioii /1Q'7^^ 
olee & Muinaii (ly/uj 


?5U Mriz 


3 /± i J 


Karuppusamy, Stappers 
&Serylak(2011) 


i iU - i(5U iVlrlz 


JU± J 
34+5 


Taylor &Hugenin (1971) 


147 MHz 


27.8 


Weltevrede et al. (2006) 


326 MHz 


30 ± 8 


Herfindal, Rankin (2007) 


327 MHz 


32±2 


Nowakowski (1996) 


430 MHz 


32 


Backer (1973) 


606 MHz 


37 


Weltevrede et al. (2007) 


1.4 GHz 


33±5 


Our Analysis (HRFS) 


8.35 GHz 


32±2 



The LRFS does not, however, resolve the issue of whether the ap- 
parent spectral feature is indeed the true one or is an alias of another 
feature outside of the unaliased frequency range. The HRFS, on the 
other hand, can help resolve such ambiguity to a large extent, and 
therefore we believe the estimate of the frequency from the HRFS 
to be more reliable. Moreover, the HRFS analysis provides addi- 
tional clues about the nature of modulation. 

Previous studies reported an LFF at fn ^ 0.033P^' (see ref- 
erences in Table 3). In our HRFS spectra (see Fig. 3), we also 
clearly see this feature as the most prominent one. There is also 
another significant feature at frequency 0.967/'[', which is cer- 
tainly the counterpart of the LFF (0.033/'[') (see also Fig. Al in 
iHerfindal & Rankiril l2007). We thus believe that the true frequency 
f^ is indeed about 0.033 P\^, and the roughly similar intensities 
(within uncertainties) of the two spectral features are to be inter- 
preted as suggesting significant amplitude modulation rather than 
phase modulation (which would be the characteristic of drifting 
subpulses). In any case, the prominence of the LFF in compari- 
son with its harmonics indicates a rather systematic (non-random) 
and possibly smooth variation of intensity across the pattern of sub- 
bearnS; 

iGil. Melikidze & Zhand ( |2006|. GMZ06 henceforth) claimed 
that, in the case of PSR Bl 133+16, the LFF reported in the litera- 
ture with the period of (33 ± 2))Pi corresponds to the carousel cir- 
culation time P^. This claim was based on the ir Partially Screened 
Gap model tail, Melikidze & Gepperll l2003h . which predicts a 
precise correlation between P^ and X-ray thermal luminosity L, 
from the hot polar cap heated by back flow bombardments - 
see Eq.3 in GMZ06, and for more details see iGil et al.l jZOOSi) . 
Based on the aveiilable X-ray Chand ra observations of B 1133+ 16 
jKargaltsev, Pavlov & Gar mire 2006). GMZ06 found the predicted 
value of Pi, = (27+2)^*1, well within the expected range (see Ta- 
ble 3fl 




' Based on the same ideas, iGil et alj j2007l) interpreted successfully the 
period (20 ± l)Pi (foun d in the fluctuation spectrum of PSR B0656+14 by 



IWeltevrede et alj|2006h as the carousel circulation time P4 in this pulsar. 



Pulse Number 



Figure 4. An example of the longitude-resolved modulation pattern (in 
grey-scale), obtained from 512 pulses folded at the putative carousel ro- 
tation period of 28.44Pi , is shown in the main panel, with the associated 
integrated pulse profile on the left. Bottom plot shows the corresponding 
intensity-modulation pattern averaged over longitudes, clearly showing an 
evidence for periodic nulls. The dashed line represents the 4cr detection 
threshold. 



This interpretation of P4 was later confirmed independently 
by HR07, although the authors admitted that they strongly doubted 
it at first. However, they surprisingly found this periodicity with 
a value of ^4 = 28. 44^1 in data with very short duration nulls, 
which they called pseudo-nulls. Such nulls can occur when the line 
of sight misses the emission between adjacent subbeams (empty 
sightline traverses through the beamlet pattern that occur occasion- 
ally but periodically). For the model which GMZ06 proposed, the 
value of P4 should not depend on the observing frequency. Our 
measurements for all available segments of data seem to support 
that claim, as the LFF-based periods we measured are very much 
consistent (within their errors) with those found from pseudo-nulls 
at lower observing frequencies. This means that the carousel circu- 
lation time stays constant over the observing frequencies ranging 
from 1 10 MHz to 8.35 GHz, almost two orders of magnitude! 

The periodic pseudo-nulls found by HR07 are among the best 
manifestations of the underlying carousel model responsible for the 
drifting subpulse phenomenon in pulsars. We have also checked our 
data for their presence, although given its erratic nature (Fig. 1) we 
were not very optimistic at first. However, to our surprise we man- 
aged to detect them in some selected, brighter pulse sequences. A 
good example from our results can be seen in Fig. 4, which shows 
the pattern of the modulation of the pulsar signal, folded with a pe- 
riod of P4 = 28.44/"! found by HR07. Each folding was performed 
over 512 pulses, i.e. about 18 putative periods P4. In four of such 
foldings (out of 16), we found a clear evidence for periodic pseudo- 
nulls in the form of apparent "null zones" in the folded modulation 
pattern. One of our best examples, as shown in Fig. 4, can be com- 



This was an especially attractive interpretation as the thermal X-ray emis- 
sion in this case obviously originated at the hot polar cap heated by back- 
flow bombardment (DeLuca et al. 2005). But what is even more important 
for us, the radio single pulses of PSR B0656+14 looked almost exactly as 
eiTatic as those of PSR Bl 133+16 at 8.35 GHz presented and analysed in 
this paper. 
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pared with Fig. 7 in HR07. One could argue, that given the erratic 
nature of our data, such nulls can appear in our folded profiles by 
pure coincidence, as we do not detect the majority of the incoming 
pulses. But given the fact that we detect about 20% of pulses, the 
chance of appearance of such coincidental periodic genuine nulls 
(as opposed to periodic pseudo-nuUs) is only 0.8'* = 1.8% (to 
minimize such chance occurrence we were forced on folding of 
a very long pulse sequences). Not to mention the fact that the fold- 
ing procedure itself helps to detect pulses that individually remain 
just below our detection threshold. Indeed, during folding of 18 se- 
quences, we improve our sensitivity by a factor of more than four 
( VTS). Therefore, we believe that the detection of periodic pseudo- 
nulls in our 8.35 GHz data confirms directly the findings of HR07 
at the much lower frequency of 327 MHz. At the same time, it ul- 
timately proves the carousel model for PSR Bl 133-1-16, that is 23 
subbeams completing one full circulation in exactly 28.44/'i . It is 
worth emphasizing that we detected our pseudo-nulls at exactly this 
period, and the LFF in our best data sequences is also close to this 
value (Tables 1 and 2). 

Summarizing, PSR B1133-(-16 seems to be a case of a pulsar 
that shows an organized behaviour of its individual pulses only for 
some limited-duration periods of time, and otherwise appears to be 
highly chaotic. We still do not know what triggers the switch from 
chaos to order, and we do not undrestand what is happening with 
the pulsar during the chaotic periods - this definitely requires fur- 
ther study. However, wc seem to have enough evidence to claim 
that, when the order appears, the pulsar radiation can be explained 
by means of the rotating carousel model. Moreover, our results pre- 
sented in this paper show that the most fundamental feature of the 
carousel model, that is the carousel circulation time (P4) prevails 
despite various difficulties nature (in form of details and nuances 
of pulsar emission) puts in our way. Whether this is some aliasing- 
smeared drifting subpulse phenomenon observed at lower frequen- 
cies, or low flux density at high frequencies preventing the majority 
of individual pulses from reaching the detection threshold, the low 
frequency feature is often visible. Additionaly, one of the most rep- 
resentative features of the carousel model is the phenomenon of 
periodic pseudo-nulls, and it is very important to note that we de- 
tected them at exactly the value of 28.44Pi reported by HR07 at 
lower frequencies. 
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